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Abstract 

The pressure for the Imperfect (Mean Field) Boson gas can be derived in several ways. 
The aim of the present note is to provide a new method based on the Approximating 
Hamiltonian argument which is extremely simple and very general. 
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1 Setup 



Consider a system of identical bosons of mass m enclosed in a cube A C l d of volume V 
centered at the origin. Let Eq < E± < E% < E% < . . . be the eigenvalues of h A := -A/2m 
on A with some boundary conditions and let {0^} with I — 0, 1, 2, 3, ... be the corresponding 
eigenfunctions. Let ai := a(<f>i) and a\ := a*((f)f) be the boson annihilation and creation 
operators on the Fock space JF A , satisfying [a;, a*,] = dij/. Let Ta be the Hamiltonian of 
the free Bose gas, that is Ta = YuHo-^i^ wnere N t = a*ai. Let N A = YlHo-^i ^ e the 
operator corresponding to the number of particles in A. The Hamiltonian of the interacting 
gas known as the Imperfect or Mean Field Boson Gas is 

H A = T A + ±N* (1.1) 
where a is a positive coupling constant, (see e.g. [T]). 

Let p := \im A ^d Eq. Let po(p) and po(p) be the grand-canonical pressure and density 
respectively for the free Bose gas at chemical potential p < po, that is: 

VM = - f ln(l - e-^)F(d V ) and p (p) = f -^JL-^Aj), (1.2) 

F being the integrated density of states of h A in the limit A j R d . Let p c := lim^-^ po (/•*)■ 
The grand-canonical pressure of the Mean Field Boson Model with Hamiltonian (jl.lj) is 

Pa(ji) = -jy In trace exp{-(3(H A - pN A )} (1.3) 

and we put 

p(fjL) = lim p A (p). (1.4) 

ATM d 

Proposition. The pressure in the thermodynamic limit p(p) exists and is given by 

( i 

-ap 2 (p) + p (p - ap(p)) if p < p c ; 

(1.5) 

(p- Po) 2 | / \ . 
2~ + Po(/io) if P > Pc 

where p c = po + ap c and p(p) zs t/ie unique solution of the equation p = po(p — ap). 

This result, for special boundary conditions, can be proved in at least three ways 0111111], 
see also El El • The aim of the present note is to provide yet another but extremely simple 
and very general way of proving this result. Our method covers also the case of attractive 
boundary conditions. The proof is based on the Approximating Hamiltonian technique, see 

e.g. m 

We shall need the following auxiliary operators for p G R and 77 e C, 

H A (rj) =H A + W (r,a* + n*a ) , (1.6) 

and 

H A (p, n)=T A + apN - \p 2 V + W (na* + V *a ) , (1.7) 



P(p) 



so that 
Let 



H A ( V ) - H A (p, rj) = ^(N A - Vpf. (1.8) 
p A (r],p) = In trace exp{-(3(H A (r]) - pN A )} (1.9) 

and 

p A (p,r],p) = — In trace exp{-/3(H A (p,T]) - pN A )}. (1.10) 
We can write H A (p, 77) — pN A in the form 

H A (p,r])-pN A = J2 e t^i^ + ^y(va* + V*ao) ~ ^ap 2 V (1.11) 
z=o 

where 

ef(p,p) :=^ A -/i + ap. (1.12) 
For convergence in (jl.K)j) must have ap — p + Eq > 0. 

2 The Proof 

The proof of the Proposition consists of four straightforward lemmas. The idea is to show 
that, for 77 7^ 0, the pressure p A (r],p) in the limit coincides with p A (p,r],p) minimized with 
respect to p (Lemmas Q and OJ) and that in turn this minimization can be performed after 
the thermodynamic limit (Lemma EJ). The final step is to switch off the source 77 to obtain 
the limiting pressure p(p) (Lemma||. 

Lemma 1 For a given rj, there is a compact subset of ((p — po)/a, 00), independent of A, 
such that for A sufficiently large the infimum of p A (p,rj, p) with respect to p is attained in 
this set. 



Proof: From ()1.10|) we can see that 



1 I I 2 1 

p A (p, 77, p) = -— \2 {In (1 - exp(-/3ef )) } + M- + -ap 2 . (2.1) 



and 



-—{p,ri,p) = -— > — 7-TT\2 +a P- 2 - 2 

dp V exp(/3ef ) - 1 (e£) 2 

Suppose 77 7^ 0. If ap < p — p + 5, with 5 > 0, then 

<9p A/ . a|r/| 2 a\ V \ 2 r , , 

^(j>,rj,fi)<-T^ + ap<-j^+fi-fM) + 6<0 (2.3) 

if 5 is sufficiently small and A large enough. On the other hand for ap > p — p + a, with 

a > 0, 

dp A , . a ^ 1 4a|77| 2 . . 4a|77| 2 

— - p,7pp > -— > — — s — +ap> -a p c -l — +ap> 

ap 1/ ^ exp(p(£// 1 — (p — &)) — 1 ar ar 

(2.4) 



for p and A sufficiently large. Therefore there exists K < oo, independent of A such that 
the infimum of p\(p, r), p) with respect to p is attained in [(p — p + 5) /a, K]. 

□ 

Suppose the infimum of p A (p,r],p) with respect to p that is attained at p A (r],p), which is 
not a priori unique. 

Lemma 2 If rj ^ 

lim Ph{r],p) = lim p A (p A ,r),p). (2.5) 

ATM d AtM d 

Proof: Since by Lemma [TJ pA is a interior point of ((// — po)/a, oo), it satisfies 

^(PA,^) = 0. (2.6) 

Now 

'-Na\ <9p A , 1 1 |?7| 



1 1 |r/| 2 



Comparing this last equation (|2.7|) with ()2.2|) we see that pA satisfies the equation 
Pa = — ( N A) HA (p A ,ri)- B y Bogoliubov's convexity inequality, see e.g. [Bj, 



< p A (p A ,r],p) -Pa(t),p) < —A A (r)). (2. 



where 

A A (r ? ) = a<(iVA-Vp-A) 2 ^ A( _ Ai??) . (2.9) 
We want to obtain an estimate for A A (r)) in terms of V. Since 



Mzl = ^r exp(^) 
aV ' V ^ (eX p (/ feA) - l) 2 + 



2 



2J77 



we use e x /(e x - 1) < 2(1 + l/x) for x > and ef (p A , p) > E$ + ap A - p. > E$ - p Q + 8 > \5 
to obtain 

52„. /o \ 1 00 1 A UI2 



g PA 
<9p 2 

By (J22I) and (J2IED 



(2 \ 1 1 4 I I ^ 

^'i^ NW-D ^W (2 - 10) 

1 °° 1 



and 

\v\ 2 



r A\2 



< Pa- (2.12) 



Therefore 

32 



' h/ ' A (p A ,r ? ,p)<2(^+/?)p A . (2.13) 



dp 2 wnr,-~y s 



Thus since by Lemma HJ p\ < K, 



Jim ^A A (r?) = Jim ^^(pa,V,») = 0. (2.14) 

V^oo V z V^oo V Op z 



□ 



Lemma 3 



1 \ri\ 2 

lim p A (r], p) = -ap 2 (r], p) + p (p - ap(r], p)) + — (2.15) 

A|M 2 ap(rp /i) - p 

\ri\ 2 

where p(rj, p) is the unique solution of the equation p = p (p — ap) + — . 

(ap - pY 

Proof: Let p(f],p) be a limit point of Pa G [(p — p Q + 6)/a, K]. p (p — ap) and p (p — ap) are 
convex in p. From (j2.1|) one can see that p\(p, T), /-0 is convex and thus as A ] R d , pa(p, r], f) 
converges uniformly in p on compact subsets of ((p — po)/a, 00) to 

1 w 2 

-ap 2 +p (p - ap) H . (2.16) 

2 ap — p 

Similarly from (12. 21) . one can see that -J— (p, V, p) + J?| is also convex and thus by the 

dp (e£) 2 

same argument, it converges uniformly in p on compact subsets of ((p — po)/a, 00). Since it 

a | 12 Q 

is also clear that A converges uniformly on the same subsets, so does — — (p, r/,p) with 
(e ) op 

limit 

a\i]\ 2 

ap - apo(p - ap) - — . (2.17) 

(ap - pY 

Since 

lim p A (r), p) = lim p A (p A , r], p) (2.18) 

ATK d ATM d 

and 

^(p A ,r ? ,p) = 0, (2.19) 

the lemma follows immediately. 

□ 



Lemma 4 



lim Pa(p) = lim lim p\(w,p). (2.20) 



Proof: By Bogoliubov's convexity inequality one has 

1 (oq + a* )) HA{lM) I <p A (p) -PA(v,fj) < "4=1 ( a o + a o)) HA ( v ,n) I- i 2 - 21 ) 



v ~ ~ JV 



< p A (p) -p A (rpp) < S (a* ) HM | < ^ (a* a )i AM < ^ (iV A )l A(7?)M) - (2.22) 



Now, by convexity with respect to p: 



iVA\ 

v / 



< Pa(v, P + 1) - Pa(»7, p) < Pa(»7, P + 1). (2.23) 



Since 



H A (7})-(ti + l)N A = T A + ^-(^ + l)N A + (VVf ] + a* )(VVr ] + ao)-a* ao-V\r ] \ 2 

> T A + ^-^ + 2)N A -V\ V \ 2 

= t a -(, -i)n a + ^(n a - v ^-^ + ^ 2 



2V V a 
V|p 



2 F(p - /i + 3) 2 



2a 

> r A -( Mo -l)iV A -VH 2 - ^~ 2 /i a + 3)2 , (2.24) 



we have the estimate 



p A fa, p + 1) < p (p - 1) + 1 + l^l 2 + (/i ^ +3)2 (2.25) 
for large A. Thus the righthand side of ()2.22|) tends to zero as r] tends to zero. 

□ 

Finally as p — > 0, p(rpp) tends to (p — po)/o if /x > /x c and to the unique solution of 
p = p (p — fl p) if p < Pc, proving the Proposition. 

□ 



3 Remarks 



Note that the proof in [2] is based on the equivalence between the canonical and the grand- 
canonical ensembles produced by the mean field interaction, aN A /2V, and the use of explicit 
information about the grand-canonical free Bose gas. 

In the proof presented in [3| the authors employ the device of shifting the spectrum to change 
the density of states removing the phase transition. 

Thirdly, in [I] a probabilistic approach based on the theory of Large Deviations is used to 
treat the pressure for this model, amongst other things. 

The essential feature of our proof is the addition of sources to make the critical density infinite 
so as to eliminate phase transitions. This allows us to linearize the mean field interaction 



and to control the particle number fluctuations. The limiting pressure is calculated in this 
regime and then the phase transition is restored by the removal of the sources. 

Finally we remark that if we replace the interaction term in (jl.ljl by the operator V<p(N\/V) 
where is a non-negative smooth convex function, the arguments used here go through with 
minor modifications to obtain the limiting pressure. 
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